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ABSTRACT 


A physical description of the destabilizing 
mechanism responsible for certain pilot induced 
lateral-directional instabilities is included 
herein. In addition, the theoretical and experi- 
mental effort propcrting to support the physical 
arguments is included. 

The basis for this study was a pilot induced 
instability known to exist for the X-15 research 
aircraft-pilot combination for certain conditions 
of flight. Accordingly, a major portion of the 
effort was concentrated upon the analysis of a 
specific known instability. Nevertheless, it is 
the opinion of the authors that the resultant 


conclusions ere quite general. 





SECTION I 
INTRCDUCTION 


The investigation upon which this report is based 
was concerned with the study of a type of pilot induced 
instability of the lateral-directional short period mode. 

In the more usual case, an instability of Gira 
nature can be explained as a function of the character- 
istic response time of a human being, in which case the 
problem is quite simply that the human controller can- 
not react to the error signals with sufficient speed. 
Quite obviously, a problem of this nature is closely 
related to the frequency of the mode of motion. 

Some recent studies, however, have indicated the 
possibility that certain pilot induced - or closed 
loop - instabilities can occur which do not appear to 
Seestrictly a function of frequency. “% Study. 
Askenas and McRuer, (Ref. 1) based upon root locus 
analyses of a pilot-airplane combination, showed the 
possibility of an instability of this type. This 
study further indicated certain combinations of the 
aerodynamic stability derivatives believed to be 
unfavorable. 

More recently, an analysis was maae by Taylor 
(Ref. 2) as a result of simulator studies of the X-15 


research aircraft. These simulator studies predicted 
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certain flight condition boundries within which the 
X-15 closed loop system would exhibit an unstable 
lateral directional short period mode. (Dutch roll 
mode). That this region did in fact exist was verified 
by flight tests. 

Both of these studies relied heavily on the root 
locus technique for an explanation of how the magnitude 
of pilot initiated control movements could destabilize 
an otherwise stable system. Consequently, althougn the 
studies showed that the pilot could destebilize the 
system, the actual destabilizing mechanism was not 
explained. This study was initiated in an effort to 
provide this explanation. 

This study was conducted in three separate, yet 
complementary, phases. The theoretical phase consisted 
of an examination of the basic equations of motion and 
the significance of certain parameters upon the resultant 
dynamic behavior or a given pilot-airplane system. Root 
locus analyses formed a major portion of this phase. 

What has been called the computer-simulator phase 
was essentially a study of the same parameters investi- 
gated theoretically, but with the aid of an electronic 
analog computer. In addition, the analog computer was 
used as a fixed base simulator for the investigation of 


an actual closed-loop system. 
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The flight phase consisted of flight tests in- 
corporating the variable stability NAvion aircraft 
described herein. This phase was primarily initiated 
in an effort to provide as rezlistic a simulation as 
possible of an unstable pilot-airplane combination. 

In general, it may be stated that the results of 
all three phases contributed significantly to the con- 


clusions which resulted from this investigation. 
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The variable stability flight simulator used in 
this investigation was a North American NAvion airplane. 
To provide the variable stability capability, a modified 
Minneapolis-Honeywell autopilot, USAF Type E-12, had 
been installed. Modificstions to the autopilot consist- 
ed of the addition of supplenentary rate gyros and a 
Sideslip transducer. Signals from tne transducers, 
proportional to the measured quantities, were introduced 
as input signals to the autopilot, resulting in control 
deflections proportional to the measured quantity. This 
made the system capable of modifying the airplane stabili- 
iy derivatives. Data collection and recording were 
accomplished by means of an ASCOP pulse width, frequency 
modulated telemet=ring system. An analog computer was 
used in the theoretical development of the problem, and, 
in conjunction with a control stick-operated potenti- 


ometer, provided fixed-base simulation. 


iio t ATRELAKE 

The NAvion is an all metal, low wing, single-engined 
monoplane. The engine, a Continental E-185, drives a 
variable pitch Hartzel propeller and is rated at 185 
horsepower — maximum continuous at sea level at 2300 rpm. 


The control surfaces are of conventional design. A 





== 
streamlined static balance is fixed at tne outboard end 
of each aileron. The trim tabs for the aileron and 
rudder are of the fixed-bend type, and tne elevator 
trim tab is adjustable from the cockpit. The physical 
characteristics of the NAvion are listed in Table I. 
A photograph of the test airplane is shown in Fig. l. 
The airplane wis equipped for this investigation with a 


three-axis variable feedback autopilot. 


AUTOPILOT 

The basic autopilot is capable of altering the 
stability of the airplane in the pitch, roll and yaw 
axes through three channels of operation. For this 
investigation, the pitch channel was rendered inoper- 
ative, with the exception of an unmodified elevator 
control. Each channel is composed of an AC series 
Summing network which combines signal voltages from the 
sensing transducers, the electrical flight controls and 
the knob actuated trim controller. The difference be- 
tween this combined signal voltage and a voltage feed- 
back signal from the servo-drum position transducer, an 
error signal, is applied to an amplifier, phase dis- 
criminated and transferred into servo drum rotation by 
means of dual power relays. The servo drum rotation 


moves the control system of the aircraft. The motion of 
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the control surface continues until error signal voltage 
at the input to the amplifier is reduced to a value 
below the threshold of the system. This velue is about 
tite, millivolts. 

The basic autopilot had been modified by deleting 
some of the original provisions and by supplying ad- 
ditional feedback loops. The directional coupler feed- 
back loop, the flight control stick trim provisiena ge a 
the coordinated turm provisions had been eliminated. 
oideslip. angle and roll rate feedback loops and a cross- 
control signal proportional to aileron control stick 
deflection had been added to the yaw channel. Sideslip 
angle and yaw rate feedback loops had been added to the 
roll channel. The manually adjustable ratio potentiometers 
were relocated on a gain control panel placed between 
the pilots! seats. An electronic safety device was 
installed to prevent "hard-over" signals by automatically 
disengaging the autopilot if an error Me ana exceeding 
approximately 3 volts should instartaneously appear at 
the amplifier input. 

Quantities to be telemetered were taken from the 
appropriate terminals on the main terminal board to a 15 
channel filter-conversion unit which converted the voltage 
to DC. A gain control for eacn channel permitted the 
DC voltage to be adjusted to that required by the tele- 


metering system (0-5 volts). From the filter-conversion 
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unit, the DC voltages representing the measured quanti- 
ties were taken to an auxiliary terminal strip and 
connected to the telemeter transmitter unit. 

The electrical controls are located in the right 
seat of the NAvion, and the standard manual controls in 
the left seat. Reference 3 described a spring centering 
device and fluid dashpots which had been incorporated 
in the electric stick. These had provided objectionable 
feel characteristics, and, according to the recommendations 
in Reference 3, were removed prior to the flight test 


phase of this investigation. 


INSTRUMENTATION 

The transducers in the feedback loops of the modified 
autopilot provided signals to the autopilot amplifier 
through the AC summing networks and provided the signals 
required to measure the aircraft motions. The flight 
condition data (airspeed, altitude, and outside air 
temperature) were obtained from standard aircraft 
instruments. Errors in these instruments were assumed 
negligible. The physical characteristics of the com- 
ponents of the instrumentation system are listed in 
table II. 

Two of the modified rate gyros were located on the 
autopilot chassis on the main equipment table and two 
were mounted on the floor of the aircraft in the equip- 


ment compartment. 
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The use of a single sideslip vane was possible 
because it contained two separate potentiometers. The 
gideslip vane was installed on a boon extending four 
feet ahead of the wing leading edge at the wing tip. 

It was assumed thet this distance was adequate to mini- 
mize measurement errors due to the wing pressure field. 
The position of the sideslip vane was such that the yaw- 
ing rate of the aircreft about the vertical body axis 
would affect the sideslip angle measurement slightly, 
but this effect was neglected. 

The telemeter transmitter unit consists of a 
rotary switch sampling each of 43 input channels plus 
two synchronizing channels at the rate of 20 times per 
second. The sampled data are converted to pulse width 
form by a keyer unit and transmitted as a UHF frequency 
modulated signal to the telemeter ground station. The 
filter-conversion unit limited the date measurements to 
14 quantities plus a full scale reference voltage. By 
jumper wiring on the auxiliary telenetering terminal 
strip it was possible to sample a given quantity more 
than once for each revolution of the switch. Specifi- 


cations of the telemeter transmitter unit are listed in 


Mable Il. 


GROUND STATION 
The ground station data collection equipment con- 


sisted of an ASCOP M series PW/FM ground station, an 
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Ampex Model 5309C dual track tape recorder and associated 
graphical recorders. 

The ASCOP ground station received the modulated 
Signals from the airborne unit, demodulated and decoded 
the signals to provide for each of the 43 channels of 
information, and provided a continuous voltage output 
to the recorders representing the measured flight data. 

The tape recorder had provisions for simultaneous 
recording of telemeter and voice transmissions. This 
permitted the recording of description of runs, pilot's 
comments and ground station operator's comments along 
with the telemeter data. 

Sanborn four channel pen recorders, Model 154~-100B, 


were used in the graphical recording phases. 


ANALOG COMPUTER 

The analog computer used in this investigation was 
a Goodyear Aircraft Corporation Model L3 (GEDA) linear 
electronic differential analyzer. Twenty-four auto- 
matically stabilized DC computing amplifiers were avail- 
able with open-loop gain greater than 5 x 10!, and of 
negligible drift. The computer incorporated an auto- 
matic error indicator and had a guaranteed accuracy of 
one percent. Provisions were available for accurately 
setting computer board potentiometers by the use of a 


special calibration potentiometer and null indicator. 
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Attempts to "fly" the computer were conducted by 
means of the simulator shown in Figure 2. This consisted 
of & Gonventional control stick which, when Gefilecez eae 
varied the voltage at the Sq terminals of the analog 
computer circuit by changing the setting of a potenti- 


ometer. 
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Computer - Simulator 


Analog 
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SBCTION III 
PROCEDURE f 


A. In-Flight Simulation 


In-flight simulation of the pilot induced lateral- 
directional oscillation was accomplished with the 
Variable stability NAvion. It was expected that results 
obtained in this manner would provide a more suitable 
basis for a physical explenation of the phenomenon than 
results obtained from a fixed base simulator, since the 
full range of visual and kinesthetic motion cues would 
be present. 

Reference 2 points out uncontrollable combinations 
of Mach number and angle of attack encountered with the — 
X-15 research aircraft witn lateral stability augment- 
ation off and tne pilot attempting to control bank angle 
in a normal manner with ailerons alone. Using these 
conditions as a guide, values were obtained for the 
dimensional stability derivatives to be reproduced in 
the NAvion. These are listed in Table III. The calcu- 
lation procedure is shown in Appendix A. These deriva- 
tives were dimensionalized for a dynamic pressure of 
w> 3 los. /ft.*, corresponding to one ¢ flizhti.)) ie. 
choice of a relatively low dynamic pressure was made in 
order that the capabilities of the autopilot/airplane 


system not be exceeded; it was anticipated that the 
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frequency of the Dutch roll oscillations would be 
reduced, while the character of the pilot-induced 
instability would be retained. Simulation at higher 
dynamic pressure was confined to the fixed base simulator. 

Dimensional derivatives for the NAvion were 
obtained for 6500 ft. standard density altitude and 
120 miles per hour true airspeed (Ref. 4). Autopilot 
feedback gain constants were then determined which would 
alter these derivatives, making them equal to those 
chosen as typical for the X-15. This procedure is shown 
in Appendix A. 

Feedback gain constants of the autopilot could be 
altered by adjustment of potentiometer settings in the 
cockpit of the NAvion. In order to determine gain 
potentiometer settings which would provide the required 
constants it was necessary to calibrate the systen.: 

Bach autopilot rate gyro was calibrated by removing 
it from the airplane and remounting it on an electrically 
driven turntable which could be rotated at a controllable 
rate. In this manner control surface deflections pro- 
duced by given angular velocities could be established. 
The deflections were measured by a propeller protractor 
for various gain potentiometer settings. 

The sideslip vane was calibrated by moving it through 


a specified angle and simultaneously measuring the control 
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surface deflections produced at different gain potenti- 
ometer settings. Electric stick calibration Weseep 721. 
€éd in the same manner, a propeller protractors7-i4. 
used to measure stick deflection. 

The results of the calibration are shown in Figure 
3, as curves of Feedback Gain Constant versus Gain 
Potentiometer Faceplate Setting. The numbering of the 
constants corresponds to the nunber of the potentiometer 
by means of which each may be adjusted. The cockpit 
panel arrangement is shovm in Figure 4. 

In order that a quantitative interpretation of the 
welenetered motion histories of the aircraft couldybe 
made, it was also necessary to calibrate the system 
for telemetry. This was performed in exactly the same 
manner as the calibration described above, with the 
exception that the measured output quantity in this 
case was voltage at the telemetry transmitter package, 
rather than control surface deflection. It was decided 
that the following quantities would provide sufficient 
information for the analysis to be performed: roll rate, 
yaw rate, stick deflection and sideslip angle. fhe 
results of the telemetry calibration are shown in 
Pigure 5 as curves of tne foregoing quantities versus 


Telemetry Input Voltage. 








a ae 
4 


5 
i 


AIO> po 
OY 





os 


1U0IS O7) ulD 











(Geer 


NAVION N9OISo6 
COCKPIT POTENTIOMETER PAiiiaam 


Ssideslip into 
yaw channel 


yaw rate 


into yaw chan. 


Sideslip into 
roll channel 


yaw rate into 
roll channel 


elevator into 
pitch channel 


rudder into 
yaw channel 


note: feedbacks 


roll rate into 
yaw channel 


aileron into 
yaw channel 


roll rate into 
roll channel 


, | 
aileron into 
Ka (Or) Kia | roll channel 


© 
3 


utilized in this investigation 


are enclosed by dashed line 















eer ee = a 
a te eee ee + meen — ee es 


| 
: 


ee ee 


FIG. 5 
TELEMETRY CALIBRATION 
= ah 3 es a I962 _ 


" 


“| 


Oo 
A 
— 
O 
a 
= 
| 
a 
eas 
| 
| 
: 
= 
| 
ia, 
- rad /se 
es Soemase 


? = 









electric stic 
deflecTrion — 





mele 

Longitudinal feedbacks, (pitch angle, pitch rate, 
velocity and angle of attack) were not employed in this 
investigation. The elevator gain potentiometer was 
adjusted to a value which gave a satisfactory elevator 
response. 

Certain deficiencies existed in the equipment which 
made it difficult to achieve an exact simulation. Among 
these were control cable stretch, structural elasticity 
of the aircraft, inaccurate resolution due to finite 
autopilot “peck sizes,” very small potentiometer settings 
and the difficulty of accurately reproducing 4 potenti 
ometer setting. As a result of these deficiencies it 
was found that the simulation procedure outlined in this 
section provided only a first approximation to the desired 
aircraft response. It was necessary, therefore, to 
further adjust the cockpit gain potentiometer settings 
in order to achieve the desired results. 

This arbitrary correction to the calculated potenti- 
ometer settings would necessitate matching of aircraft 
responses to analog computer responses if the objective 
was exact simulation. A reasonably good qualitative 
simulation of such motions can be readily achieved, 
however, and was considered to be sufficient for the 


purpose of this investigation. 
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BP, Analog Computer Study 

Three basic tools were used to accomplish the 

computer phase of this investifation. These were: 

1. the GDA Analog Computer, 

2. the Sanborn Recorder, and 

dD. the Control Stick Simulator, 
The computer and recorder were used to do the open Loop 
studies of the system, while all three units were used 
to complete the closed loop work. 

The equations mechanized in the computer are listed 
and discussed in detseil later in tnis report. 

When open loop response to step inputs by the 
computer matched calculated values as to time to damp 
to 1/2 amplitude and period of oscillation, attempts 
were made to "fly" the computer. These closed loop 
trials were accomplished by using the control stick 
Slmaulator in the same manner as the control stick in an 
aircraft. The time histories of bank angle, roll rate, 
Sideslip angle, and yaw rate on the recorder were used 
es references. The control stick simulator produged 
input voltages to the computer analogous to aileron 
inputs to the basic airplane. 

Three separate configurations of the basic equations 


of motion were investigated. These were: 
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A. ®& = O and low dynamic pressure (153 psf). 
Bis Oe 


10° and low dynamic pressure (153 psf). 


Cage es 


10° and high dynamie pressure (1000 psf). 


The wiring diagrams for the computer study are 


shown in Pigures 6 and 7. 
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DISCUSSION AND RESULTS 


mhneoretvical Development: 


During the past several years, considerable emphasis 
has been placed on the determination of what parameters 
can best be used to measure the lateral-directional 
handling qualities of an airplane. As a result, a consider- 
able body of literature has come into being and certain 
airplane response characteristics have come to be re- 
garded as unfavorable or undesirable from the standpoint 
of control by a human pilot. (Ref. 1). Primarily as a 
result of the human pilot being involved as a part of the 
feedback control loop, investigations of this nature have | 
mot been easy to interpret either physically or mathe= 
matically. 

This particular investigation has been concerned 
with a rather special segment of the lateral-directional 
handling qualities problem ~ specifically the case where 
the inclusion of a pilot in the control loop has resulted 
in an unstable Dutch roll oscillation, although the open 
loop Dutch roll response of the aircraft was known to be 
stable. In short, this investigation has been concerned 
with the specific problem of a pilot induced Dutch roll 
instability. 

The problem of a pilot-airplane combination of the 


nature described above could have been regarded as general 
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and an overall solution sought. However, it was felt 
by the authors that a more concentrated effort could be 
made if the problem were approached from the opposite 
point of view - namely a detailed theoretical and 
experimental investigation using, as a guide, a pilot-— 
airplane combination for which the results were known. 
The reasoning followed was that a detailed study of a 
particular situation might then yield a theory which 
could be investigated in a more general sense. 

According to Ref. 2, the X-15 research airplane- 
pilot combination was found to result in an unstable 
lateral-directional oscillation for certain flight 
conditions if the stability augmentation system was made 
inoperative. For tnis reason, the X-15 airplane-pilot 
combination was the closed Loop system chosen for the 
detailed study. Specifically, the instability associated 
with a flight Mach number of 3 and a ten degree angle of 
attack was investigated. 

The fact that the ‘above configuration was the one 
chosen for study placed several important restrictions 
upon the problem and thereby upon the resultant conclusions. 
First, altnough the above configuration was handled 
throughout this study as though the aircraft had been in 
steady flight, in reality the X-15 encountered the above 


stated flight condition only in transitory flighty. 
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Accordingly, the aerodynamic stability derivatives were 
continuously changing, with the result being that the 
pilots of the X-15 were advised to refrain from using 
the rudder as a control. Therefore, this entire investi- 
gation was conducted assuming the ailerons to be the only 
type of lateral-directional control available , although 
the effect of using rudder is discussed in the flight 
phase section of this discussion. A second condition 
implied by the problem as stated is that the frequency of 
the oscillation was not a factor. That this must be so 
can be reasoned by noting that the frequency is more a 
function of the dynamic pressure than of Mach number, 
while the instability has been described as being a 
function of Mach number and angle of attack. Thus, the 
solution sought was not one dependent upon the inability 
of a human being to respond with sufficient speed to an 
error signal. In effect, a suitable explanation to the 
problem encountered by the X-15 would have to be valid 
even if the pilot were assumed to respond without a time 
lage 

The stability derivatives required for a detailed 
study of the above condition were obtained from Ref. 5 
and are listed in Table IV. These derivatives lead to 
autopilot potentiometer gain constants listed in Table V. 


It is noteworthy that only two of the stability derivatives 
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were Significantly different from what one might ordi- 
narily expect. Namely, both Chg, (aileron yaw) and Ca, 
(dihedral effect) had a positive sign. Thus, these two 
derivatives immediately became suspect and it was expected 
that a theoretical analysis of the problem would show 
either one or both of these derivatives to be involved 
in an explanation of the instability. In this rather 
special case of a single control, pilot~induced insta- 
bility, the unusual sign of the aileron yaw became of 
Particular interest. 

A proposed theoretical explanation to the problem 
encountered by the X-15 aircraft was published by Taylor 
(Ref. 2). In his work, Taylor showed, by root locus 
techniques, that the Dutch roll mode of an aircraft with 
the characteristics exhibited by the X-15 could be made 
unstable if a pilot gain proportional to roll angle and 
roll rate was included in the feedback loop. According 
to Taylor, the ratio of “/wy could be used to define 
the region where an instability might occur - a ratio 
greater than one being indicative of an unfavorable 
condition. The human transfer function used by Taylor 
in his analysis was determined empirically and found to 
be K(1 + .57s). Thus a pilot response proportional to 
both roll angle and roll rate was assumed. The analysis 


followed by Taylor was quite believable for two reasons. 
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First, the transfer function was seen to be independent of 
a pilot time lag which would seem to isolate any effect of 
frequency from the basic control problem. Secondly, the 
emalysis as carried forth by Taylor resulted 1m Viegeom. 
clusion that the basic cause of the problem was the positive 
sense of the dihedral effect CC ga). What Taylor's analysis 
did not do was present a physical explaration of the phe- 
nomenon, 

In an effort to assure that the problem had been 
interpreted correctly by the authors, the initial theo- 
retical effort was directed toward matching the results 
of Taylor's investigation. Thus the human transfer 
function proposed by Taylor was used and a root locus 
study of the effect of varying pilot gain undertaken. 


The set of equations used was as follows: 
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It is worthy of note that the above set 15 [oneag aa. 
Cipal axes coordinate system. This particular coordinate 
System was chosen for this study inasmuch as the stability 
derivatives as obtained from Ref. 5 had been calculated 
with respect to this system of axes. 

The result of the root locus study can best be 


explained by considering separately the effect of varying 
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the pilot gain, K, and the trim angle of attack, Wp. 
The fact that &, was a significant parameter actually 
was determined during the analog computer phase of this 
program and is discussed in detail in connection with 
that part of the program. Here the fact that O&,is import- 
ant will merely be stated and then substantiated by the 
root locus technique. 

Assuming the pilot gain to be the variable of interest, 
the determinant of the coefficients of (1) was expanded 


with the following root locus equation resulting: 


(2) (A+Ar)(A= As) CA® + 2EyWnyA+ wry’), uCldaK. 0 
(AE + 2EpWngpA + Wng*)(\ + 57 AY) ur 
The derivation of the above equation is included in 
Appendix B. 
Figure 8a was the type of root locus which resulted 


from the above equation when & owas assumed to be negligible. 


This root locus showed the inclusion of a proportional 
pilot gain to be a stebilizing fector at zero angle of 
attack, provided the stability derivatives could be 
assumed constant. Furthermore, a study of the parameters 
which determine the value of Wn¢g, (Sq. 5 in Appendix B), 
reveals that an increase in either of the two suspect 
stability derivatives (Cg, and Cig) would actually de- 
crease the value of Wng and thus tend to stabilize the 
closed loop system. This is a result of the signs of 
both derivatives being reversed from what might be termed 


normal. 
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Ba. Root Locus for a&.=O , O<K< 00 


The next step was to include the correct value of My 





(10 degrees). From equations (1) it can be seen that the 
inclusion of @, has no effect upon the denominator term 
of the root locus equation but does alter the character- 
istic equation somewhat. The added terms which are the 


result of including G, are: 


oe 
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It can be seen from equation (3) that with a positive 
Pig and all other derivatives of the usual sign, the 
effect of including & is to decrease the A and Ne 
terms of the characteristic equation. This had little 


effect upon the roll and spiral modes but did result in 
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both the damping ratio (Ey ) and the undamped natural 
frequency (Wn ) of the Dutch roll mode being appreci- 
ably decreased. The root locus which resulted is depicted 
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Fig 8b. Root Locus for a, 10°, OC K<c 


This is similar to what Taylor obtained. Furthermore, 





it can be seen from equation (3) that the destabilizing 
effect of Xo can be attributed to the positive sign of 
the dihedral effect. (Cg). 

Having shown the importance of the trim angle of 
attack, it was deduced that a root locus using trim 
angle of attack (@,) as the varying quantity might be 


informative. This study was carried out by assuming the 
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stability derivatives to be independent of the trim 
angle of attack. Thus only the orientation of the air- 
craft with respect to its flight path was assumed to be 
variable. The root locus equation for this investigation 


had the following form: 


(4) (A+ drllA+ As + CE yWny +Wny*) 1 7 
MA+ KK) a 


with the denominator terms found from equation (3) and with 


mee Ce ee = 
K ry dK 5T 9, (CagCn, ae CnpCe,) 


Figure 8c shows the root locus for O£&,Soofor no 





pilot gain while figure 8d shows the same result for an 
arbitrary pilot gain. Calculations for the conditions 
investigated showed the system with pilot gain to go 
unstable for a trim angle of attack on the order of 
twenty five percent of whet was required for instability 
without the pilot in the loop. Conversely it was noted 
that for a trim angle of attack of ten degrees, the open 
loop system was so slightly stable that a small pilot 
gain could be expected to be sufficient to cause the 
closed loop instability. From equation (4) it can be 
seen that a negative Chie would change the sign of K’* 

and in turn change the angle condition for the root locus 
diagram. For this case, the effect of increasing Ae 


would not be destabilizing. 
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several interesting observations resulted from a 
study of the root loci briefly discussed above. First, 
it was observed that if the signs of both the dihedral 
effect C Coa) and the aileron yaw (Cnoae) had been 
reversed, there would have been virtually no change in 
the closed loop response at low angles of attack. How- 
ever, a study of ene cofactor of the ((,) term showed 
that a beneficial effect would have been realized as 
the trim angle was increased. Secondly, the observation 
was made that for a positive dihedral effect, the pro- 
verse sense of the aileron yaw definitely tended to 
improve the closed loop stability while an adverse 
aileron yaw would have had the opposite effect. That 
is, an adverse aileron yaw would have placed the complex 
zero (Fig. 8b) even further above the complex pole. 
This was believed to be particularly significant inas- 
much as with pilot time lag removed as a significant 
contributor to the instability, one might reason that 
the only way that the pilot could be contributory as a 
destabilizing factor would be as a means for introducing 
an unfavorable forcing function while actually responding 
correctly to the error signal. In the case of a single 
control, this meant an unfavorable input into either the 
rolling or yawing moment equations by means of an aileron 
deflection. For the particular case investigated, the 


root locus plots showed a proverse aileron yaw (CnSa ) 
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to be favorable while an increase in the trim angle (,) 
was seen to have a destabilizing effect. Thus it was 
reasoned that the scorns physical effect was that an 
increase in trim angle of attack opposed the proverse 
aileron yaw effect until, at some angle of attack, the 
sense of the sideslip resultant from an aileron deflection 
was probably reversed. This might be explained by noting 
that a simple uncoupled rolling motion about the body 
axis of an aircraft will result,as the roll angle increases, 
in the conversion of angle of attack into sideslip angle. 
For example, a positive roll would tend to produce a 
positive sideslip. This is the opposite of the effect 
produced by a proverse aileron yaw. Thus it is reason- 
able to assume that at some point the effect of angle of 
attack and of proverse aileron yaw might be equal. At 
this point and beyond, an increase in the pilot gain 
would cause the pilot to introduce a destabilizing force. 
Looked at in the light of Taylor's analysis, this point 
where the aileron yawing moment due to aileron deflection 
reversed direction is probably coincident with the point 
where W$/Wy becomes greater than unity. 

In summary, the theoretical effect agreed quite well 
with Taylor's analysis although it is felt by the authors 
that the approach involving the variation of trim angle 
of attack provided a basis for the physical reasoning 


explained above. However, it can be seen that even if 
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the physical effect of increasing the trim angle is as 
stated, the mechanism by which the instability occurs 
cannot be explained from the theorv alone. Fortunately, 
a study of the time histories of both the analog computer 
Simulator and the variable stability NAvion did provide 
an insigcht into the mechanism involved and did tend to 


agree with the hypothesis suggested by the theory. 


Analog Computer - Simulator Phase: 


As mentioned previously, the approach used in this 
investigation was to make a detailed study of an airplane- 
pilot combination which results had shown to be unstable. 
For the computer phase of the investigation two variations, 
along with the basic condition, were examined. The basic 
condition was the X-~15 research airplane flying at 
M = 3.0 with an angle of attack of ten degrees. The two 
meriations were: to keep all other conditions ttemean— 
but decrease the trim angle of.attack to zero, and to 
keep all other conditions the same while increasing the 
dynamic pressure from 153 to 1,000 pounds per square foot. 


All three of these conditions will be discussed below. 





For this condition, the basic equations of motion, 
using the non-dimensional X-15 coefficients, reduce to 


the following set of computer equations: 
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Both step aileron inputs and "“attempts-to-fly" the 
computer simulator showed this condition to be completely 
manageable without‘ the slightest tendency to go unstable. 

This fact, in conjunction with the other studies 
outlined herein, confirmed the significance of the trim 
angle of attack (&,) to the study of the closed loop 


instability. 





Condition If: M = 


This condition resulted from an analysis of the 
mesuits of Condition I. By re-evaluation of the basie 
equations, using @,= 10 degrees, the following set of 


equations was obtained: 


100% = -3.358 - 100 + 1.06T¢ + 17.45 ¢ 
IOP = GOB 41.7250 - 2.126+ 83.9dc 
200 = 55% ~ .9024 -.0902f + 9.08da 


These equations were mechanized into the computer and 
resulted in the predicted pilot-induced instability. 
A step aileron input and an "attempt-to-fly" the 


computer are shown in Figures 9 and 10. From the step 
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input response there were several important things noted. 
First, it was seen that the motion was lightly damped. 
Second, there was considerable excitation of the Dutch 
roll mode from a coupling of the rolling and yawing modes 
of motion. Third, altiough the sign of Cn$8q was posi- 
tive - leading one to believe the yaw, and consequently 
the sideslip, due to aileron should have been proverse - 
the total effect of the parameters affecting the side- 
Slip was just the opposite. That is, the sideslip angle 
response described by these equations, due toa coo 
aileron input, was as might be expected from an adverse 
yawing moment. The fourth and lest item to te noted was 
the phase ansle difference between the angle of roll and 
the angle of sideslip. It could be seen from the time 
histories, thet the two motions were almost exactly out 
Ot phase. The signiricance of the phasing of theses aase 
while the Dutch roll motion was in progress, as well as 
the sign of the sideslip resulting from an aileron input, 
will be discussed more fully later. 

Using the control stick simulator, "attempts-to-fly"” 
the simulator resulted in oscillations of increasing 
amplitude if the "pilot" responded in a normal manner 
to bank angle and roll rate. Using the stick to make 
aileron inputs in an effort to keep the wings level 
resulted in an increase of the sideslip angle , due to 


the ailsron deflection. Thus, sideslip grew, fed the 
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roll excursions, and the resultant situation was one of 
a closed loop instability. By releasing the stick, the 
oscillations damped out, although the roll angle grew 
indefinitely as a result of the unstable spiral mode. 

An "attempt-to-fly" the computer by using the a 
technique described by Petersen in Reference 6 was 
completely successful. The oscillations were allowed to 
build up using normal pilot techniques, and then, at 
some point, ailerons were used primarily to stop the 
sideslipping, and Fe eends to keep the wings level. 

As the 3 trace came through zero, the stick controller 
would give the stick a sharp pulse in the direction the 
trace was moving. In the airplane this would correspond 
=O a sharp right aileron pulse as the nose of the alrpiane 
went through zero degrees sideslip heading toward the 
right. In this manner, the motion could be damped out 
within a few cycles. It did not appear feasible to 
maintain steady state conditions using this technique, 
but oscillations resulting from normal pilot inputs could 
be rapidly damped out in this manner. Roll angle Conwro 
was maintained by varying the magnitude of the pulse and 
the stick position between pulses. 

During the study of this condition, the following 
observations were noted: 

1. The oscillatory motion was lightly damped, 

2. There was considerable coupling between 


the roll and sideslip motions, 
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3. The sideslip generated by an aileron deflection 
was adverse even thoughCnSq was positive. 

4. During Dutch roll oscillations the angles of 
bank and sideslip remained approximately 180 
degrees out of phase. 

5. Normal attempts-to-fly were unsuccessful, and 

oe The technique could be used to damp out any 
Dutch roll oscillations that developed. 





®o= 10 degrees 


This condition was investigated primarily to see if 
the characteristics of the motion were affected by changes 
in dynamic pressure. The basic equations of motion modi- 


fied to increase the dynamic pressure were written as: 


20 8 = -4,.5223 - 20 + 1.44 3 3.49o 
2D = 116.88 42.32 - 2.84% +108.4é5q 
4p = T1.2B6-1.212%-.1212 + 11.12 dq 


These equations were mechanized in the computer and 
both step aileron inputs and "attempts-to-fly" were 
observed. The basic character of the motion did not 
change. All the comments recorded under Condition II 
are valid for the higher dynamic pressure., The only 
significant change in the motion was the frequency of 
oscillation of the Dutch roll mode. Previously the 
motion had a period of about five seconds. At this in- 


creased dynamic pressure the period was reduced to 
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approximately two seconds. This did not alier the diffi- 
culty of attempting to fly the configuration. Nor did 
the increased frequency invalidate the technique of 
Sonera . 

This condition can be summed up simply by saying it 
was identical to Condition II except things happened a 
little faster. However, the motion was not so fast as 
to make tne controller's task too difficult, nor so fast 
that the character of the motion and subsequent responses 
to the motion were changed. 

One result of the above observations was that the 
frequency of the motion (within the limits investigated) 
was seen to be an insignificant contributor to the in- 
stability. As long as the motion proceeded at a pace 
which the controller could follow, and as long as the 
normal control motions were made, the basic instability 
occurred. Thus the frequency of the motion was elimi- 
nated as an important parameter insofar as a study of 
the basic destabilizing mechanism was concerned. This 
confirmed the hypothesis outlined in the theoretical 


mortion of the discussion. 


Blight Phase: 
The flight phase of this investigation involved the 


use of the variable stability NAvion. The physical 


characteristics of this airplane, its simulation equip- 
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ment and the accompanying instrumentation are described 
elsewhere in this report. 

The initial flight phase was based upon an in- 
stability encountered by the X-15 researcn aircraft in 
a steady flight (lift = weight) condition at Mach 3 and 
at ten degrees angle of attack. The dimensional X-15 
derivatives used in calculating the potentiometer 
settings are listed in Table III. A direct method for 
Simulating the X-15 side force equation was not available, 
and the dynamics introduced by the autopilot equipment 
were not known by the authors, hence it was necessary to 
correct the computed autopilot potentiometer gain 
settings in flight until a desired condition was achieved. 

Flight simulation was directed toward achieving the 
desired character of pilot induced instability, rather 
than toward an exact simulation of the dynamic response 
of the X-15. On this premise, and on the basis of 
computer studies, it was decided that the most signifi- 
cant characteristics to be simulated were the following: 
1) a lightly damped Dutch roll mode 2) a roll rate 
response to an aileron function which contained a signifi- 
cant Dutch roll effect. 

Cockpit potentiometer settings were adjusted until 
the two above cnaracteristics were obtained. It was 
found that the period of the Dutch roll mode obtained 


in this manner was approximately three seconds, and the 
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steady state roll rate of the aircraft was slightly de- 
creased. Thus Duten roll period is of the same order 
of magnitude as the periods of the Dutch roll modes 
studied on the computer. In addition, computer studies 
indicated that the frequency of the Dutch roll mode was 
unimportant insofar as the fundamental problem was 
eoncerned. 

Undoubtedly the most important characteristic was 
Simulated, namely, the closed-Loop pilot-airplane system 
exhibited an unstable Dutch roll oscillation in spite of 
the fact that the open loop response was licshtly damped. 
Thus, as in the case of the X-15, the attempts of a 
pilot to improve upon a marginal condition had a de- 
Stabilizing effect. Significantly, this unstable result 
could be avoided if the pilot used both ailerons and 
rudder for control. Just why the authors feel this 
point to be significant will be discussed later in this 
meport. 

Time histories of the various angles and rates as 
recorded on a typical flight are included herein. In- 
asmuch as the mode of interest was the Dutch roll, the 
authors felt the two most significant time histories to 
be those of the sideslip angle (6) and the roll rate 


(~). Accordingly, the time histories of these two 


items were studied in detail and will be the items 


discussed herein. 
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It can be seen from the time histories oz the air— 
plane response to a step aileron input (Fig. 11), that 
moe roll rate did contain a considerable amount Gievares 
moll, and further, that the Dutcn roll motion wasiizns 
damped. Also, it can be seen that the proverse sense of 
the aileron yew resulted in the initial sideslip generated 
by an aileron input being of opposite sign to the roll 
rate response. This was as might be expected for an air- 
plane with appreciable proverse aileron yaw and only a 
Slight amount of yaw due to roll. However, this initial 
Sideslip response to an aileron input was the opposite 
of what was observed to happen on the analog computer 
Emautator. With the benefit of hindsieht, it canvse 
stated that this difference was prohably due to the fact 
that the effect of the trim angle of attack (Q,) had 
not been considered. 

From the time histories of the attempts to fly the 
resultant configuration, (Fig. 12) it can be seen that 
attempts by the pilot to control the oscillation resulted 
in excursions of increasing magnitude for both the side- 
Slip angle and the roll rate. Releasing the controls, 
(Fig. 13) on the other hand, resulted in decreasing 
amplitudes of motion. However, once large amplitudes of 
motion had been obtained, what appeared to be a limit 
cycling of the autopilot prevented 2 complete damping of 


the open loop motion. An important piece of information 





acre al 
RESPONSE TO 
STEP AILERON 
INPUT 


Navion Variable 


Stability Aircraft 
Basic X-I5 Settings 





— +20 
—+20 
0 y deg 
c. | sec 
) ee 2O 
Fa —-I5 
pa mt | YW O ds, deg. 
| —+|5 


, deg 





—-|5 
0 4 8 |2 time, sec 





Das | out} 


bap' of 





— Sit 


| a 


FOUN ALGELS eqdMmHA UOIADKY 
ALITIGVLSNI GSONGNI - Lod 
cl Olas 





FIG .13 
OPEN LOOP 
DUTCH ROLL 
CONVERGENCE 


Navion Variable 


Stability Aircrat+ 


O 1O 20 +ime , sec. 








-+|5 





=O , deg 


--15 
~ piloted — stick free 





aoe 
that can be obtained from these tine histories is the 
phase relationship of the sideslip angle ( ) and the 
roll angle (Q@-) - the latter being assumed to lag the 
roll rate by a quarter of a cycle. In particular, a 
can be seen that the sideslip angle and the roll angle 
are almost exactly in phase. This was also in contrast 
to the phasing obtained by the electronic analog simulator. 
Thus it was noted that both the analog simulator and the 
NAvion simulator exhibited a pilot induced oscillation 
in spite of the fact that the character of the motion was 
virtually opposite for the two cases. That the sideslip 
resulting from a step aileron input was also different 
for the two systems suggested the possibility that, in 
spite of the apparent difference, the same destabilizing 
mechanism was involved in each case. Accordingly, the 
aileron yaw was reversed from the proverse to the adverse 
sense and the tests repeated. Time histories (Fig. 14) 
show that although the response of sideslip angle to an 
aileron input was reversed, the percent of Dutch roll 
in roll was maintained at a high level and the Dutch 
roll remained lightly damped. However, this configuration 
could be controlled by the pilot. This is not to say 
that the resultant configuration was a good - or even 
acceptable - airplane from a handling gquelities stand- 
point. The damping of the Dutch roll mode was so slight 
that unpleasant oscillations developed rather easily - 


particularly if rapid control reverssls were attempted. 
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However, it is significant that the pilot was able to 
control the aircraft, and, in fact, damp the oscillations 
when they occurred using only normal aileron control 
movements. This seemed to confirm the hypothesis enee 
the destabilizing mechanism for the X-15 pilot-aircraft 
combination (as manifested by the an2rlog simulator) was 
of a similar character to what was obtained with the 
PAvion. In particular, it could be seen (that wine 
stability resulted when the pilot, in an attempt to 
control roll angle, succeeded in reducing the roll angle 
only at the expense of increasing the sideslip angle. 
During this flight phase, flight of the paca 
(closed loop unstable) configuration using both ailerons 
and the rudder to maintain lateral-directional control 
was tried. The results showed that by using the rudder, 
the sideslip angle could be controlled and the aircraft 
handled comparatively easily. Rapid reversals presented 
a problem as it was difficult to keep the aircraft in 
balanced flight throughout the maneuver and a rather 
large Dutch roll oscillation often resulted. However, 
by using rudders, the oscillations could be damped within 
a few cycles, (Fig. 15). This also tended to confirm the 
importance of the sideslip angle to the oscillation. 
Aerodynamic derivatives other than the aileron yaw 
were varied and a pilot opinion of these variations 


obtained. However, even though many of these derivatives 
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were found to have an important effect uron the handling 
qualities, none were found to be so directly involved in 
the actual destabilizing action of the pilot as the 
aileron yaw, but rather were indirect contributions to 

the problem. As such, the effects of these other stability 
derivatives were probably of importence in establishing 

the conditions believed to be necessary for a problem 

of this type to exist. These conditions are further 
discussed in the following section. 

The reason for the difference in phasing of roll 
and sideslip between the computer and NAvion is believed 
to lie in the inability to directiy simulate the side 
force equation with the iiAvion. An attempt was made to 
adjust the NAvion stability derivatives in such a way as 
to allow for the effect of trim angle of attack (Wp). 
Flight test of this eltered configuration, however, 
indicated that the phasing had not heen changed. 

The fact that the NAvion did not yield the sane 
phasing of sideslip and roll as was obtained with the 
analog computer was not considered by the authors to be 
a shortcoming of the flight simulation. It wes viet 
rather, that this difference in phasing helpeqwa. 
emphasize the more feneral nature of the physical 
explanation of the instability mechanism proposed here- 


in - that is, since the same character of pilot induced 
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instability was manifested in both of the two different 
types of systems investigated, the problem is not limited 
to one particular aircraft. Furthermore, since the 
analysis of the mechanism of instability is consistent 
with both of the two types of motion phasins, the appli- 


cations of the analysis also become less limited. 


Sumnary of Results: 


Quite obviously the different phases of this program 
have resulted in a multi-sided approach to the basic 
problem of a pilot~induced Dutch roll instability. In 
general, it can be stated that each phase of the investi- 
gation contributed something unique. Thus, in this section 
of the report, these separate and distinct contributions 
will be summarized in an effort to show how each con- 
tributed toward the formulation of what is believed to 
be a correct physical explanation of the destabilizing 
phenomenon. 

From the root locus investig2tions discussed pre- 
viously, it was determined thet, with a positive dihedral 
effect, an increase in the trim angle of attack (MQ) 
had a destabilizing effect upon the characteristic open 
loop response of the X-15. Further, a study of the 
equations of motion showed that the effect of increasing 
the trim amgle of attack was to increase the coupling 


of the sideslipping and rolling motions of the aircraft. 
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In addition, a study of the co-factor of the trim angle 
of attack (Q@,) term, in conjunction with the root locus 
Shown in Figure 8c, indicated that the phasing of the 
motions was of such a nature that destabilizing moments 
were generated. Even more significantly, the root locus 
analysis showed that, under certain conditions, the normal 
inputs of a pilot attempting to keep the wings level 
had a destabilizing effect upon the Dutch roll mode. 
It was then reasoned that inasmuch as pilot response 
time had not been considered in the analysis, the de- 
stabilizing effect of the human controller could not 
have been the result of erroneously applied corrections 
to unwanted deviations in roll mgle, but rather had to 
be due to the indirect effects of the resultant aileron 
deflections. This indirect effect was seen to be the 
conversion of argle of attack into sideslip. Thus it 
was concluded that the attempt of a pilot to control 
roll angle resulted in the generation of a sideslip 
angle. Further, the phasing between the rolling and 
sideslipping motions was of such a nature that one 
motion tended to reinforce the other. 

The results obtained as a consequence of the analog 
computer-simulator phase of this program complemented 
the theoretical investigation and, in addition, provided 
an explanation of the actual destabilizing mechanism. 
The computer studies confirmed the importance of the trim 


angle of attack to the oscillatory behavior of the system. 
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Subsequently the root locus investigation was under- 
taken and the results discussed above obtained. 

Study of the time histories obtained during the 
computer phase of the investigation showed that the 
Dutch roll oscillation, for the equations represent- 
ative of the X-15, was characterized by angles of roll 
and sideslip which were of opposite sign. FPurther, 
the open loop response to a step aileron input showed 
that an input producing a positive roll resulted in an 
initial positive sideslip (Fig. 9). This was not what 
was expected, inasmuch as the aileron yaw was proverse 
and the yaw due to roll was small. However, as has been 
previously explained, the conversion of angle of attack 
into sideslip as the result of a rolling motion could 
produce an effect of this type. 

By combining the two effects noted above - which 
would be the natural result of a pilot attempting to 
control a Dutch roll motion with ailerons - it could be 
seen that an attempt to control the angle of bank ina 
normal manner would result in an input tending to in- 
crease the angle of sideslip. With the high degree of 
coupling present between sideslip and roll, as evidenced 
by the considerable Dutch roll contribution to the roll 
rate time history, an unstable behavior resulted. This 
hypothesis was checked by artificially increasing tle 
proverse aileron yaw to such an extent that the initial 


sideslip generated by a step aileron input was changed. 
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Thet is, an input which resulted in a positive roll 
yielded a negative sideslip. Figures 16 and 17 show 
tne time histories resultant from tnis configuration. 
With the increased proverse aileron yaw, the configuration 
could be controlled in spite of the fact that the Dutch 
roll mode remained lightly damped and the percentage of 
Dutch roll in roll remained aprreciable. The lightly 
damped Dutch roll mode did result in easily induced 
oscillations but, unlike the basic configuration, pilot 
attempts to control roll angle did not destabilize the 
system. In fact, as might be expected, it was found that 
attempts to control the roll angle in a normal manner 
also helped to reduce the excursions of the sideslip 
angle. 

A by-product of the computer study was the confirm- 

ation that the frequency of the oscillation was not a 
Significant contributor to the basic instability. Quite 
obviously, this does not mean that frequency can always 
be dismissed as insignificant, but only that an instability 
of the type studied could exist even at frequencies 
completely within the control limits of a human being. 
In this phase, it was found that motions with a period 
as high as five seconds were completely unmanageable 
as long as a normal pilot technique was used. 

The flight phase of this investigation yielded 


results which essentially complemented those obtained 
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ED 
from the fixed base simulator. This was in spite of the 
fact that the craracter of the oscillatory motion obtained 
with the NAvion variable stability airplane was consider- 
ably different from what was observed with the fixed 
base simulator. 

As outlined previously, the connection seemed to be 
that, in either case, an attempt by a pilot to control 
the roll angle by manipulation of the ailerons resulted 
in an unfavorable generation of sideslip. Inasmuch as 
the rolling and sideslipping motions were highly coupled, 
the increased sideslip generated additional rolling 
moments which in turn required even greater control 
movements. The result was a coupled rolling, sideslipping 
motion of ever increasing amplitude. That the oscillations 
could be controlled - or even damped - by judicious use of 
the rudder, substantiated the validity of tne foregoing 
explanation of the instability mechanisn. 

The result of the combined Studies Was )tl2teeene. 
combinations of roll angle, sideslip angle, and side- 
Slip resultant from an aileron deflection, were found 
to be unfavorable. In particular, it was noted that 
the pilot-induced unstable behavior occurred only when 
the pilot, in an effort to return to a wings level 
condition, could not reduce the bank angle without 
increasing the sideslip angle. Thus an “effective” 


proverse aileron yaw was found to have a destabilizing 
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effect if the Dutch roll oscillation was of such a nature 
that the bank angle and the sideslip angle were in phase. 
For an “effective” adverse yaw, the opposite phasing was 
unfavorable. Multiple exposure photographs depicting 
the two motions described above are included herein. 
(Pigs. 18 and 19). The term "effective" as used refers 
to the sense of the initial sideslip produced by an 
aileron deflection. This initial sideslip depends on 
the combined effects of Cnéq »>Cnp and Q&,, rather than the 
actual sign of the yawing moment produced. Thus it might 
be termed an integrated effect. 

In addition to the above stated phase relationships, 
certain conditions were necessary before the instability 
could occur. First, it was necessary that the lateral- 
directional motions be highly coupled. Physically this 
was indicated by the large oscillations which occurredyin 
the time histories of roll rate response to a step aileron 
input. Secondly, it was necessary that the Dutch roll 
mode be lightly damped. No formal proof of this last 
statement has been offered herein, but it seems quite 
obvious that a highly damped mode would virtually never 
demand pilot attention. Practically speaking, the pilot 
would probably never be aware of its existence. Accord- 
ingly, one could hardly expect the pilot to ever initiate 
the destabilizing mechanism described above. From a root 


locus standpoint, an increase in the Dutch roll damping 
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would shift the Dutch roll poles and zeroes (Pigs. 8a 
and 8b) to the left. This would obviously have a 
stabilizing effect. 

From the discussion of the destabilizing mechanisn, 
it is apprrent that the phasing between the bank angle 
(p) and the sideslip angle (() is an important indicat- 
or of how a closed loop instability might develop. In 
particular, this phase relationship determines the sense 
of the effective aileron yaw that can be regarded as 
unfavorable. Accordingly, a study was undert2ken in an 
effort to ascertain what parameters determined this phase 
relationship. This study is outlined in Appendix C. 


Briefly, the phase angle between P and was found to 
“Ndq 
dQ 


This same combination of derivatives has been found to be 


be heavily dependent upon the sign of the term:-2uCygesea 


important in the determination of the 2, parameter. 

Further analysis led to the conclusion that the 
phasing of the Dutch roll in roll could be expected to 
bear some relationship to the character of the roll rate 
response. Figures 9 and 16 showed that a distinct 
difference existed between the characteristic time 
constants of the roll mode for the favorable and un- 
favorable situations. In particular, the time difference 
was seen to be approximately equal to one half of the 


Dutch roll period. This effect is discussed in detail 
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in Appendix D. Briefly, a shift of 180 degrees (Dutch 


roll in roll) was seen to coincide with a shift through 


Weg 
Oy 
In sunaary, for the type of pilot induced instability 


unity of the magnitude of the parameter. 

outlined herein, tne arguments presented appear to be a 
physical way of describing the phenomenon depicted by 

the root locus studies of Taylor. The theoretical 

effort conducted in conjunction with this investigation 
seems to show that the onset of an unfavorable combination 
of phase relationships is coincident with tne parameter 
geccribed as “kK, becoming greater than unity. 

Although only a limited number of cases of pilot- 
induced instabilities were studied in detail in con- 
junction with this investig-tion, the authors have seen 
no reason to suspect that the destabilizing mechanism 
as outlined herein cannot be regarded in a broader sense. 
Thus it is the opinion of the authors that an explination 
of this type might well be valid for other instabilities 


known to be a function of pilot gain. 
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CONCLUSIONS 


As a result of this investigation, the following 
conelusions were drawn: 
1. The NAvion variable stability aircraft can be 
used for a general study of pilot-induced 


lateral-directional instabilities. 


2. A pilot induced inst2bility of the type investi- 
gated can be successfully demonstrated ,and 
analyzed with a fixed base simulator incorpo- 


rating an electronic analog computer. 


?- For this type of pilot induced instabilii7ge. 
the Duten roll mode to occur, certain Condreien 
were found to be required. These are listed in 
qualitative fashion as: 

a. A lightly damped Dutch roll transient; 
b. Highly coupled rolling and sideslipping 
motions - as manifested by a large per- 


centage of Dutch roll in “rows 


4. If the conditions listed in item (3) were prevalent, 
certain phase relationships between the roll angle, 
Sideslip angle, and the effective aileron yaw were 
found to result in an unstable closed-loop system. 
In particular, if the Dutch roll oscillation snowed 


the roll angle and the sideslip angle to be of the 
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Same sign, an unstable system resulted if the 
effective aileron yaw was proverse. For the 
opposite phasing, an effective adverse aileron 
yaw resulted in an unstable closed loop system. 
For the type of pilot induced instapvia 
studied, the freauency of the oscillation was 


found to be an unimportant p* rameter. 


Even a marginally stable system (open loop) 

could be managed quite easily if a direct means 
for controlling the sideslip angle were provided, 
Accordingly, a system which was closed loop un- 
stable with only aileron control was found to be 
manageable if the rudder was used. Also, the (A 


technique could be used, although greater care 


was required. 


It was further concluded that this study provided 
a physical interpretation of the type of pilot- 
induced instability predictable by the method 
described by Taylor (Ref. 2). 
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TABLE I 
PHYSICAL CHARACTsRISTICS OF THE NAvion 


WING DATA 


Total Area (including ailerons, flaps 


and 19.87 ft. covered by fuselage) 183.54 £t.* 
Span 33, 39 ete 
Aspect Ratio 6305 
Taper Ratio 0.54 
Dihedral Angle 7.50 deg. 
Reet Chord [-2C ste 
Mean Aerodynamic Chord 62.3550 


Incidence Angle 


Root 2.0 deg. 

Tip -1.0 deg. 
Sweepback of Leading Edge 3.0 deg, 
Twist 

Geometric 3.0 deg. 


Airfoil Section 
Root NACA 4415R 
Jualye, NACA 6410R 


Flaps, 40 deg., Plain 


AIL®3RON DATA (For one aileron) 
Span 61.99 in. 


Deflection 30 dés. up; 20 deg .aamae 
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TABLE I (continued ) 


AILERON DATA (continued) 
Control Wheel Throw 
Aerodynamic Balance Frise-type nose 
Static Balance (Outboard end of 
each aileron) Streamlined weight 
Trim Tab (Right aileron) Fixed Bend Tab 
Ratio of Aileron Chord to 
Wing Chord 0.284 


HORIZONTAL TAIL DATA 


Total Area (Including 2.37 £t.° 

covered by fuselage) 43.05 £t.* 

Span 15.1 7 eee 

Aspect Ratio 4.02 

MAC Brey ir 

Airfoil Sections, Root and Tip NACA 0012 

Incidence Angle -3.0 deg. 
BLEVATOR DATA 

Total Area 14.10 £t.° 

Span 73.58 ise 

Deflection 30 deg. up; 20 deg. dwn. 

Deflection Trim Tabs 

(32 in. span, 4-1/2 in. Chord) +30 deg. 

Root Chord 1.5 £86 


Tip Chord ].O fae 





TABLE I (continued) 


VERTICAL TAIL DATA 


Area (Including 2.57 ft. blanketed by 


fuselage and excluding 1.84 ft.* dorsal 


ey 


fin area) 12.93 ft.° 
Span 4.05 £8. 
Airfoil Section 
Root NACA 0013.2 Modified 
Tip NACA 0012-64 Modified 


Incidence Angle (With respect 
to FRL) 


2 deg. Nose Left 


RUDDER DATA 
Area 6.05 £t.° 
Deflection 17 deg. L3 23 deg. R. 
Trim Tab Pixed Bend Tab 


Rig (Angle with respect to 


fin center line) 


FUSELAGE DATA 
Fuselage Length Over-all 
Width, Maximum 
Depth, Maximum 


Fineness Ratio 


POWER PLAN? DATA 


Airplane is powered by one Continental 
Maximum continuous rated horsepower at 


185 at 2300 RPM. 


3 deg. right 


27.25 ft. 
4.14 ft. 
4.40 ft. 
Gae 


E-185 engine. 


sea level 





TABLE I Ccontinued ) 


POWNR PLANT DATA (continued) 


oF 


Tne propeller is a Hartzel hydro-selective propeller 


with the following characteristics: 


Activity Factor 
Diameter 


Pitch (.75R) 


MISCELLANEOUS DATA 
Weischt 
Basic 
Fuel (40 gal.) 
Pilots (2 with parachutes) 
Gross Weight 
Center of Gravity Position 
Tail Length 
Autopilot Mechanical Gear Ratios 
‘Deg. Blevator per Deg. Servo-drum 
Deg.-Aileron per Deg. Servo-drum 


Deg. Rudder per Deg. Servo-drum 


100 
86 in. 
26.5 deg. 


2129 lbs. 
240 lbs. 
_410 lbs. 
2119 its. 
29.5% MAC 
16.663 aie 


0. 9a 
0.910 
0.330 





TABLE II 
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SPECIFICATIONS OF INSTRUMENTATION COMPONENTS 


RATE GYRQS 
Manufacturer 
Model Number 
Power Input 


Power Consumption 


Rotor Speed 
Weight 


Maximum Turn Rate 


Potentiometer Resistance 


Potentiometer Excitation 


SIDESLIP VANE 
Manufacturer 
Model Number 
Potentiometer Resistance 


Potentiometer Active Angle 


Potentiometer Excitation 


Minneapolis-Honeywell 

JG 7005A-24 

115 voltsy§00 cycle ace 
32 watts (starting), 

13 watts (running) 
20,000 RPM 

1.75 lbs. 

Modified to approximately 
30°/sec. 

530 ohms 

Gyro #9 - 24 volts A.C. 
Gyros #10, #11 and #13 - 
30 volts A.C. 


Giannini 

2516 

1997 ohms 

Yaw channel - + 50° 
Roll channel- + 45° 


Yaw channel - 15 volts A.C. 
Roll channel- 30 volts A.C. 
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TABLE II (continued) 


TELEMET2R TRANSMITTER UNIT 


Manufacturer 
Model Number 
Input Signals 


Information Channels 


sampling Rate 

Accuracy 

RF Power Output 
Frequency Range 
Frequency Stability 
Primary Power Requirement 


Weicht 


ASCOP 

DT-4 

O-5 volts D.C. 

4% plus two for synchro- 
nization of ground station 
20 RPS 

+ 1% of full scale 

4 watts 

215-235 MC 

+ 0.05% 

28 volts D.C. 

19.6 lbs. 
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TABLE IIT 


Dimensional derivatives of the X-15 at Mach = 3.0, &,= 100 


Le = +9.02 No = =e008d 

ome 6 lie NO, = +6454 

Ng = 42,75 LS = +8.40 
TABLE IV 


Non dimensional stability derivatives of the X-15 


foeetech = 5.0, &%,= 10° 


Coa = +,044 Che = -1.40 

Cy, = +280 CnSa = +.052 

Co, = +-228 CiSa = Face 

Cg = +,.3515 CL = eae 

Cnp = ~-.140 ar = -1.38 
TABLE V 


Autopilot potentiometer feedback gain constants. 


4 = +.541 deg. rudder/deg. sideslip 
Ke = -~.121 deg. rudder/deg./sec roll rate 
Ke = ~,lll1 deg. rudder/deg./sec yaw rate 
Ka = +.609 deg. aileron/deg. sideslip 
Kg = -.0756 deg. aileron/deg./sec yaw rate 
Ko/Ky 3 = -.345 deg. rudder/deg. aileron 
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Calculations for the Plight Simulation: 


Non-dimensional stability derivatives chosen as 


typical for the X-15 at M = 3.0, @,= 10° are: 


CI, = + 

CQp = = 

CQ. = + 

4 = + 
"p 


For o 


O44 Ch 


2 
~280 Cnn 
228 ove 
Jie ore 


ne - g flight: 
W = 


r 2 Cy v S 


13,445 - 


1/2 (.44)(3 x 1005)°(200) = 3.36 x 107-lb-sec*/ft* 


— 
<— 


~ 140 Cn = Cu =e 


a 


oe 


(S 
O52 


O41 


Where 1005 ft/sec is the speed of sound corresponding to 


an altitude for which 0 = 3.36 x 107°1b-sec"/#4" 


oS 1/2oV 


Z 


= 1/2 (3.36x107")(9.1x10°) = 153 1b/ft* 


-. 13,445 _ ; 
TU = M/OSV = 38257075 (200)(3015)(32.174)= 206 


; _ 13,445 _ 
P= M/OSb = 3°8e5T0~> (200)(32.174)(22.36) = 2790 


2 


2 
k_° = Ix = 3,348 x 32.174 = 8.02 ft 
X MW 13,445 
2 2 
k © = I = 78,691 x 32.174 = 188.4 ft 
“ g =  =13,445 
J = 2 = = 2 8.02 = .0321 
0 
J, =2 (ke) * = 2 188.4 = .753 


9, 


(2) 
(3) 


(4) 


(5) 
(6) 


(7) 


(8) 
(9) 
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Values for dimensional derivatives were obtained as 


follows: 
La = 55 : ,044(2790 me 6) 0. (10) 
; pert Stay 
Lp = ion = -,280 ere (11) 
BG). 220. cuiaesosel 
_. 2 Seer - 4.172 (12) 





2TI, BBO. CYT. O32T) 


Nea = am BIG 3132, = + 2675 (13) 


N45 = ae aes o1l4 = ee O04 oe 14 
P2tJ, eet0.cc7 “a 
Ne oe = -.0451 (15) 








Nbe= aan = 4052(2790)_ = +4454 (16) 
— Cte - = 
daz e499 = 041 (2790) = + 8.40 (17) 
er “ C. 


The equations which express the relationship between the 
simulator airplane derivatives and the simulated deriva- 
tives are derived in Reference 4. The resultant express- 


ions for the autopilot feedback gain constants are: 


ky = 1.312 - 2.415 x 10°7Lg@ - 0.2774Ng (18) 
ke = -0.123 - 2.415 x 1077» - 0.2775 Np (19) 
ke = -0.124 - 2.415 x 10°?Ly ~ 0.2775 Nr (20) 
Kg = 0.346 + 0.0490 La + 0.0278 Na (21) 
Kg = 0.373 + 0.04901 5+ 0.0278 Np (22) 
k,, = 0.0828 + 0.0490|, + 0.0278 Nr (23) 


10 





a 
Kis 


Be O= 


~ 20.5(Néa/|da) + 0.178 (24) 


2-05(N&q/\_da) + 3.62 


Values obtained, fron the foregoing equations, for feedback 


gain constants were: 


Ke = 


+ .541 deg. 
- .l2l deg. 
- elll deg. 
+ .609 deg. 
+ .362 deg. 


rudder/deg. sideslip 
rudder/deg/see roll rate 
rudder/deg/sec yaw rate 
aileron/deg sideslip 


aileron/deg/sec roll rate 


Kyo = ~ -0756 deg. aileron/deg/sec yaw rate 


Ko/Ky 3 = - .345 deg. rudder/deg. aileron 


Cockpit potentiometer settings were obtained from the 


calibration curves for these values of feedback gain. 
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APPENDIX B 
The following is a derivation of the root locus 
equation used to study the effect of pilot gain. 
A pilot response, proportional to roll angle, of 
the form K(1 + .57d) is assumed. 
The determinant set, in non-dimensional form, is 


as follows: 


Bd e bo 
Cyp 2d a2 Craeoea O 
HCNe Cle, Sled - Jxd* BOL ae 
Cn Cnng - Jd aly MCnSq 

O O K(1 + .5Td) — | 


Expansion of the above can be symbolized in determinental 


form as: 
Cya- 2d sae Ci Tae 
Che Ca,./2 _ Gee 
BCne Cnr - Jed <REg : | 
(2) 
hee ae O 
+ K(1+57a) | eCep = Pe/2 uChsa | = O 


BCnia Cnn ~ Jd MCnéa | 
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kxpansion of the first determinant set will yield 

a quartic equation which can be divided by ed J, and then 

factored into a form from which the characteristic modes of 

transient motion can te deduced. Thus this is the character- 


istic equation and will be symbolized as: 


ID] = (A+AR(A+ ASM AR + 2 EyWny + Wny*) a 


The remaining determinant can be expanded and the 


sum represented as follows: 
ID] +K a og [pCa 2)(CR i Br Jad) + a pCnel 
~CnSq (vg 2d) ous + 2 uci 


eae, 


This can be rearranged to yield: 


ID| + K ( (\+ BPC iio fae (Chet JaCva)a + CY Ete png) 
CU, 25 (5) 


Cnda_ Cor Ic 2. 1a) 
+ d- | Lye + <4=* 
CLdaq CJ oF Clr 
The terms multiplying the human transfer function are 
seen to be a quadratic which is known to yield a complex 


pair of roots. Accordingly the above may be symbolized 
as follows: (dis replaced by) ) 


JD] + wOlbaK (1 + STAAS + 28 pug * Ong") (6) 


The root locus equation is thus found to bes: 


(A+Ar)At AslAS+ Z2SyWny A + Wn) 4 wlan 20 


i 


(E+ 28$00ng A + Cng2)(i + .57A) z (7) 
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APY ENDIX C 
Discussion of Phasing of Boll and Sideslip: 

One of the parameters discussed in conjunction with 
the instability studied herein is the phase relation of 
roll angle and sideslip angle. This relationship is the 
subject of the following discussion. 

From the basic determinant for the lateral- 
directional motion of the airplane, given here in Laplace 


operator form as: 


6 ay 6 Se 

eo —Z, . a8 CoS ; : (1) 
poe ee 

BEng Ene S Jne EDLs a 


one may form the ratio p to/o as: 











Cy_ - <9 a? O 
rr _ Cds 
pChie Che ee q 
4 BCnp one Sie - KEnoa 
i - (2) 
O =e C, + ae 
eG " O Vinee 2 
en ne 2° 
—nr _ =O 
a es 2° a 
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This may be written as: 


er: 
7) =s +As + B 
Cs +Ds+E 


Cy 


) 


where 
ih = Cnéa one Z Cy Pe. Cnr 
Cidq C.Jz Gay 
_ uCne , CyaCne _ Cnba Cy is eta) 
pa Bene , CypCnr _ CndaClr [Cye Mie 
aon . i » GEOUSa zum , . rag 
C = 2Oy. Cnda - Coogee 


Cia 
y= Crp = COR + Qo Cnr -- ; + (Cho + QLo Cir) 


et oe e0Cr,. = Cnda CLCOr 


ee ee 


Za Che 25 


For the condition of low dynamic pressure (153 1b/ft*), 
with Q&. = 10° this ratio can be written 

3° + 1.8l1s + 962.1 (4) 
ee ae ee 


See 1 S51 
ieee | Constant 


These roots are plotted on the complex plane and yield a 


ple “Ore 
oO 
Fy 


phase relation for °/3 of -173 degrees (Fig. 20). 
It appears from an examination of the above expression 
of ?/ 2 that in those cases where the Dutch roll mode is 


lightly damped, the damping terms of the numerator and 
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denominator will also be small. This means that all the 
roots will lie close to the imaginary axis. Further, 
by analyzing the plot of these roots, it is evident 
that the phasing of the o/6 motion is primarily de- 
pendent upon the relative position of the numerator 
root and the Dutch roll root. If the numerator root 
has an imaginary part greater than the imaginary part 
of the dutch roll root, the phasing between PD and 
will be approximately 180 degrees. If the Dutch roll 
imaginary part is the larger, then the phasing will be 
nearly zero degrees. 

Those parameters which affect the imaginary part 
of the Dutch roll root are difficult to isolate since 
all approximations applied to the condition under investi- 
gation yielded poor results. However, it is evident 
that a major contribution comes from BSP. This 
combination of parameters also appears in coefficient B 
in the expression for 2/6 , so its influence in shift- 
ing the two roots is equal. The only other significant 
term in coefficient B is 2 Cig Gaba » It is through 
this expression that a major shift in the numerator root 
is achieved. In the case under consideration, this one 
term islarge enough to shift the numerator root below 
the Dutch roll root, changing the phasing to zero degrees. 
This shift would occur if one of the three parameters had 


@ sign change. 





EB7 = 

This is merely additional supporting evidence for 
the argument that the phase relation between PD and {© 
and the yaw due to aileron are important parameters 
in the explanation of pilot induced lateral-directional 
instabilities. It also points out the problem which 
will arise in any airplane - pilot combination where 
Cig changes sign within the performance envelope. If 
this change of sign occurs in an area where other 


conditions are favorable to pilot-induced instabilities, 


then these instabilities may well be encountered. 
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APrENDIX D 
Discussion of Dutch Roll to Roll Mode Phasing: 

In an attempt to shed more light on the mechanism 
producing the type of instability umder consideration, 
the following analysis of Dutch roll to roll mode was 
made. Of particular interest was the phase relation 
of the two motions. 

Again the basic determinant for the lateral- 
directional motion of the airplane is given here in 


Laplace operator form. 


‘es _ oa 


Cyg- 2s -2 OL+2u5 0 
el! Gs) - 2 
uCue a Jx8"  - WCMSa 
BCne . JAS Nes - UCnda 


Forming an expression for D in terms of this 
determinant and applying a step input yields the follow- 


ing determinant ford. 


Cyg- Cs -2 O 

CIB Cir/2 = Ca 

/Cnp Cor_ Js -UCnda | 
P . Crp - 2S -2 CL 

UChe C8r/2 Chps - JS 


MCna “or - Jas nes 


This may be written as: 
Als + C) 
S- 


; “+ Bs | 
(s- A\S-dal[s- As)[S- Aa) 





where 
Apo 
oe es (Cre ea -Cya) 


— 


we a ee 


COE: CIdq - Vilage) 2 
_ (CraCne/2) + 2pCng _ Cle Cnda So 2 aC 
ames aus rE oy. 2 C0Sa a cca 


and the denominator represents the characteristic of 
motion. 
For the condition of low dynamic pressure (153 1b/f£t*) 


with Wo= 10°, this expression becomes: 


$ . 3565 (s*+ 1814s + Sol) 


(S-AS-Az\(s - Aayls- A) 


The =~ of the above are plotted on the complex 





plane (Fig. 21). It can be seen that the phase angle 
changes by 180 degrees whenever the numerator root moves 
inside the Dutch roll root. Further, a comparison of the 
roll traces (Figs. 9 and 16) shows the same result. 
Figure 9,with effective adverse yaw, shows a time delay 
of 1.5 seconds to reach the first peak on the roll trace. 
Figure 16, with the same conditions except effective 
proverse yaw, shows @ time delay of .5 seconds. This 
difference in time corresponds to one half the period of 
the Dutch roll motion. Here again, as in the PD to @ 
relation, we observe a change in phase of 180 degrees at 


the same point where We Judd becomes greater than l. 
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